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Introduction 


To understand the purpose of much of the mathematics in this unit, we need to 
look back to a group of units earlier in the course: Unit 6 (second-order 
differential equations), Unit 7 (oscillations) and Unit 8 (damping). In Unit 6 we 
studied the solution of differential equations, such as, for example, 


xX + 4x + 10x = 5cos 3t. 


In Units 7 and 8 we studied topics in mechanics where such second-order 
differential equations are used. In Unit 7 we studied oscillations, and in particular 
simple harmonic motion, the ‘fabulous perfect’ oscillations described by equations 
of the form : 


X¥ + w’*x = 0. (1) 
In Unit 8 we studied the model 
X + 2awx + w*x =0 


which takes account of the damping exhibited in real mechanical systems. The 
behaviour of such a mechanical system under an externally applied sinusoidally 
varying force was also studied. This may be modelled by an equation such as 


X + 2awx + w*x = Acos(vt + d). (2) 


All the mechanical problems studied in Units 7 and 8 concerned the motion of a 
single object in one dimension. Suppose we wish to study the motion of a 
mechanical system involving two (or more) objects, or of a single object whose 
motion cannot be described by just one position coordinate. In such cases we may 
arrive (say by an application of Newton’s laws) at a description of the motion of 
such an object involving linked differential equations in more than one variable, 
such as 


X= —3x + $y 

: Pe (3) 

y= ge a. 
In the television programme associated with this unit we will examine the motion 
in a plane of an object whose position coordinates (x,y) are governed by these 
equations. We will also see how the equations can be solved. 


The substance of this unit is mathematics, although the motivation for much of its 
study comes from mechanics. The unit covers methods for solving certain systems 
of differential equations, such as (3) above, or 


2x; + Skype x, — 4x5 + oe 


(4) 


X) = 2X, = dey, + Sink 


In mechanics we are usually concerned with systems of equations, such as (3), 
involving second derivatives (and known as second-order systems). In this unit we 
will first study the simpler case of first-order systems (such as (4)) where only first 
derivatives occur. We will consider a method of solution for homogeneous first- 
order systems in Section 1, and a method for solving inhomogeneous first-order 
systems in Section 2. In Section 3 we study the solution of homogeneous second- 
order systems of a special form, similar to (3) above. Such systems may be 
regarded as an extension of simple harmonic motion (Equation (1)) to more than 
one variable. Finally, in Section 4, we look at systems of equations generalizing 
Equation (2) (which describes forced and damped motion) to more than one 
variable. We will not study the general solution of such systems, but will seek to 
extend to systems ideas such as ‘transient’ and ‘steady-state solution’ important in 
Unit 8. We shall also see how to compute steady-state solutions for such systems. 
The steady-state solution is often all that is required in mechanics. 


You have seen how systems of linear algebraic equations may conveniently be 
written using matrices. Systems of linear differential equations may also be 
expressed in matrix form. We shall find that the ideas of eigenvalue and 
eigenvector of a matrix, introduced in the preceding unit, are of value in solving 


We discuss the modelling of 
such mechanical systems in 
Unit 24. 


These linked differential 
equations were called 
‘coupled’ differential equations 
in Unit 15. 


Although the term ‘system’ is 
usually used in this context we 
could equally well talk of sets 
of differential equations. 


The term ‘homogeneous’ is 
defined in Subsection 1.1, and 
has much the same meaning 
as for a single differential 
equation. 
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systems of differential equations. Thus in this unit you will need to bring together 
ideas from earlier units on two apparently diverse topics: differential equations 
and matrices. 


Study guide 


The sections of this unit are not of equal length. You will probably find that 
Sections 1 and 2 take longer to study than Sections 3 and 4. 


There is no audio-tape for this unit. The television programme is part of Section 3, 
and the reading specifically associated with the programme is in Subsection 3.1. It 
is not necessary to have studied Sections 1 and 2 before watching the programme; 
however, you should read the pre-programme notes in Subsection 3.1 before 
viewing the programme. You will also find it helpful (but not essential) to have 
read the introduction to this unit and Subsection 1.1 before viewing. 


When studying this unit you will need to be familiar with the following material 
from previous units: 


Unit 2, Subsection 4.2 (the integrating factor technique); this is required for 
Subsection 2.1 of this unit. 


Unit 21, particularly Sections 1 and 2 (eigenvalues and eigenvectors); this is 
required throughout this unit. 


Unit 6, Section 4, or Unit 8; this is required for Section 4 of this unit, where we 
need: 


(i) | the structure of the solution of a forced oscillation equation of the form (2) 
above; in particular, the ideas of ‘transient’ and ‘steady-state solution’; 


(ii) the calculation of steady-state solutions by the method of phasors. 


Unit 2, Section 2—reviewed in Unit 19, Section 1—(Euler’s method of numerical 
approximation); this is required for Subsection 2.2 of this unit. 


1 First-order systems 


1.1 Matrix notation and terminology 


As I mentioned in the introduction, a system of linear algebraic equations may 
also be written in matrix notation. For example 


ey = Xs + 2xu =. «5 3 -1 2| |x, 5 
xX; — 7X, + 4x3 = -1 may be written 1 —-7 4) |x| = |-1l 
Ke x= QD 0 1 1| |x3 0 


or, more concisely, 


Ax =b 
3 -l Z 
where A is the matrix |1 —7 41, x is the vector [x, X, -%,)° and bis 
0 1 1 


the vector [5 —1 oy. 


Let us now consider the possibility of writing a system of differential equations in 
matrix form. 


The system 
ae + 3X5 = X14 = 4x> + e' 
(1) 


xX; — 2X, = 4x, + Xs + sint 
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can also be written 


Fee! a Pe 


To take full advantage of matrix notation we would like to write this as 
Ax = Bx + h(t), 


2 3 1 -4 : 
where A = a B= ; ,p xis the column vector [x, x |’, and 


h(t) = [e sint |’. Here, the entries in the column vector x are functions. We have 
extended the idea of differentiation to such vectors in a natural way by writing 

x= [X, X,]". As far as vector algebra is concerned, column vectors of functions 
can be treated just like ordinary vectors. 


Example 1 
Suppose f(t) = [t? 3+t  eé)' and g(t) = [2t? ee oul 
Then 

(f + g)(t) = [327 7+t e+e]? — (addition) 


df (t 
— = [2t 1 ey. (differentiation) 
Exercise 1 


i) Fina &@ 


dt 


(ii) If x = ae”, where a is a constant vector, find x. 
[Solution on p. 32] 


, where g is defined in Example 1. 


Using such vectors of functions, we can readily write suitable systems of 
differential equations in matrix form. 


Example 2 
Write 
x, = 3X5 = Zn eoeke 
: : ; (2) 
2X, + 4X, = x, + 3x,+4+sin4t 
in matrix form. 
Solution 
In matrix form System (2) becomes 
AX = Bx + h(t) 
1 —3 2 -1 
where A = , B= eee Ey,, x2]' and 
yy 4 1 3 
n()= (0 sin4r}. 
The matrix notation need not be restricted to systems of two equations in two 
unknowns, nor to matrices with constant entries. 
Example 3 
Write the following system of equations in matrix form. 
tX; —X3=sint 
X2+X3+ tx, =0 (3) 


(sin t)X> — X1 = 2X9 
Solution 


First write the system as 
tx, —X3=sint 


Xz +X34+ tx, ae 


—X, - (sin t)x, —2x, me A}, 
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Then it is easily seen to be equivalent to 
A(t)x + B(t)x = h(t) 


t 0 0 oe 4 
where A(t)=| 0 1 fet t 0 Ol eee oh yt 
—| sin t 0 0 -—2 0 


andh(t)=[sint 0O OJ. 


However, not every system of differential equations can be written in matix form, 
in the way we have in these examples. For example the system 


7 + ee =X] + X2 
ae (4) 
X1X2 — X41 Zo X2 
could not be. 
In this unit we shall be concerned almost exclusively with systems of the form 
A, (t)x + A2(t)x = h(¢) (5) 
or of the form 
A, (t)X + A,(t)x + A3(t)x = h(t) (6) 


(where A,(t) is not the zero matrix in either case). Such systems are said to be 
linear. Equation (5) is a linear first-order system (it contains no derivative higher 
than first-order). Equation (6) is a linear second-order system. These definitions of 
‘linear’ and ‘order’ are extensions of those for a single differential equation. We 
shall usually be concerned with constant-coefficient systems, in which all the entries 
in the matrices A,(t), A,(t) and A3(t) are constants (though h(t) may depend on t). 
If h(t) = 0 in (5) or (6), then we have a homogeneous linear system. 


Exercise 2 


Which of the systems (1)-(4) in the text above are linear? Of those that are linear, say 
which (if any) are homogeneous, and which (if any) are constant-coefficient, and give the 
order in each case. 


[Solution on p. 32] 
Exercise 3 
(i) Write the system 
Xo + 2 = 3X5 —-X> 
X3 =X, +xX2+X3 
Xo =X, +X; 
in matrix form. 
(ii) Is this system constant-coefficient? Is it homogeneous? What is its order? 
[Solution on p. 32] 
In the remainder of this section we are concerned with methods of solution of 
linear, constant-coefficient, first-order systems of differential equations. In each 


case the first step in the method of solution is to write the system in normal form 
as defined below. 


A linear constant-coefficient first-order system is in normal form 
when it is written in the form 


x = Bx + k(t). 


The general constant-coefficient linear first-order system is 
A,x = 5 A>x — h(t). 


This can be written in normal form so long as the matrix A, is non-singular. For 
then we have 


% = —AT'Ax + Ay "h(?) 


which is in normal form. 


See the introduction to Unit 6. 
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Exercise 4 


Where possible, write each of the systems of differential equations below in normal form. If 
it is not possible, explain why. 


2X; —X2+ 10x, =2 X, + X23 =0 
(11) xX, = ia X + 5X1 =f (iv) X1 + ss = [X; — X2 
2X; + 2X. +x; —x,=0 X, —X2 =X, + 2x2 


[Solution on p. 32] 


1.2 The solution of constant-coefficient homogeneous systems 


We now consider a method of solution for constant-coefficient, first-order, 
homogeneous systems. 


Method of solution 


Let us consider a constant-coefficient, first-order, homogeneous system in normal 
form: 


x = Bx. (7) 


We can derive quite a simple formula for the solution of (7), based on eigenvalues 
and eigenvectors. Let us first see where the formula comes from. To do this, we 
look for solutions of (7) which have the form 


x = ae” (8) 
where a is a non-zero constant vector (and J is a constant scalar). 
If x = ae“, as in (8), then 
x = Jae” (see Exercise 1 (ii)). 
Thus (8) is a solution of (7), so long as 
jae" = Bae“ 
That is, so long as 
Aa = Ba (where a # 0). (9) 


But this equation is just the definition of the eigenvalues and eigenvectors of the 
matrix B. So we see that x = ae” is a non-zero solution of x = Bx so long as / is 
an eigenvalue of the matrix B and a is a corresponding eigenvector. 


Now the matrix B will (usually) have more than one eigenvalue, so x = Bx will 
(usually) have several solutions of the form x = ae". When we have a 
homogeneous linear system of differential equations it is quite easily proved that 
any linear combination of solutions is again a solution. So if A, and /, are 
eigenvalues of B, and a, and a, are corresponding eigenvectors, then 


x = C,a,e*" + C,a,e*# 


is also a solution of (7) (where C, and C, are arbitrary constants). 


We have found a number of solutions of'(7). The only question that remains is 
whether we can construct the general solution of (7) by taking linear combinations 
of exponentials in this way. Unfortunately the answer is ‘not always’. The matrix B 
must satisfy the conditions given in Theorem 1 below. Fortunately, many of the 
matrices which occur in practice do satisfy the conditions of the theorem. 


Theorem 1 


Let B be a square, n x n, matrix. Suppose that B has n linearly 
independent eigenvectors, say a,,42,...,a,, Corresponding to the 
(not necessarily distinct) eigenvalues /,,/5,...,A,. Then the 
general solution of the system of differential equations 


x = Bx 


x= Cae + Cia.e” +... + Cae, 


where C,,C>,...,C, are arbitrary constants. 


If a = 0 we get x = 0, which is 
a solution of (7), but not a 
very interesting one. 


You may notice a similarity to 
the general solution of a 
homogeneous second-order 
differential equation, in Unit 6. 
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Theorem 1, then, provides a procedure for finding the general solution of a 
homogeneous linear first-order system, so long as the system satisfies the 
conditions of the theorem. 


Procedure 1.2: To solve a linear, constant-coefficient, homogeneous, 
first-order system 


For this procedure to work we need to be able to express the 
system in the normal form 


x = Bx 


where B is an n x n matrix with n linearly independent 
eigenvectors. 


Express the system of equations in the normal form x = Bx. 


Find the eigenvalues 4,,A2,...,4, of B and a corresponding 
set of linearly independent eigenvectors a,,a>,...,a,. 


Write down the general solution of x = Bx in the form 
x = Cyaje"" + Cane +... + Cae 


where C,,C2,...,C,, are arbitrary constants. 


Example 4 
Find the general solution of 
X¥1=X%1+X2 | 


X = 3x1 seo 


Solution 
The system is already in normal form x = Bx where Step 1 
b 
B= 
3° =] | 
We must find the eigenvalues and eigenvectors of B. Step 2 


The eigenvalues are found by solving the characteristic equation det(B — AI) = 0. 


Now 
1-—A 1 
det(B — AI) = 
3 —1-j;j 


= (1 — A)(-1 — A) -3 
= A? — 4. 
So the eigenvalues are 2 and —2. 


To find corresponding eigenvectors we solve 


Peaks! 
= 0 
3 —1—A]| v 
putting A equal to each eigenvalue in turn. 


Case A = 2: we obtain the equations 


—u+ v=0 
3u — 3v=0 
which have solutions of the form u = k, v = k. So an eigenvector corresponding to 
the eigenvalue 2 is [1 tT. | The general eigenvector is 
[k  k]", but we only need 
Case 4 = —2: we obtain the equations one eigenvector here. 
3u+v=0 (twice) | 


sO an eigenvector corresponding to the eigenvalue —2 is [1 —3]!. 
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Since the conditions of Procedure 1.2 are satisfied (B is 2 x 2 with 2 linearly 
independent eigenvectors) we can write down the general solution: 


1 1 
x= caf fe + c,} HP 


That is — 
ye fe +e 
x, =e —3Cxe ? 
where C, and C, are arbitrary constants. 


It is sensible to check solutions to differential equations. This is readily done here. 
With x, and x, as above, we have 


ky = 20,6 = 20a 
while 
x1 + X2 = 2C,e* — 2C,e~*. 
SOX, =X, + X>,-as8 required. 
Similarly 
= 7 ie 4 OC he 
= 3x; —X>. 


Exercise 5 
Find the general solution of 


Xx; = 2X1 + 3x, 
X2 = Ix; + X2 


[Solution on p. 32] 


An example with a repeated eigenvalue 


For Procedure 1.2 to work the matrix B must have n linearly independent 
eigenvectors. It does not matter if there are less than n distinct eigenvalues 
provided there are enough linearly independent eigenvectors corresponding to each 
eigenvalue to give a total of n linearly independent eigenvectors. Let us see an 
example of this. 


Example 5 

Find the general solution of the system 
X, = 5x, + 3x3 
X2 = 3x, + 2x, + 3x; 
X3 = —6x, — 4x3. 

Solution 


The equations are in normal form, as required. We have x = Bx with 


5 
Be 3 3 3 
=—6 = =4 


We need to find the eigenvalues and eigenvectors of B. 


The eigenvalues are found by solving the characteristic equation det(B — AI) = 0. 
We have 


ee are 3 6 S23 
oi eae | a ek an eee eee Ss 3 
OS ee pee a ee 


The choice of eigenvectors is 
not unique. Other choices lead 
to different looking, but 
equivalent, forms of the 
general solution. 


See Unit 21, Subsection 1.4, 
for a discussion of 
eigenvectors corresponding to 
repeated eigenvalues 


Step 1 


Step 2 
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oe ae 3 
a 
a eo - 2-2) : 


i 
0 —-6 -4-A 


—4-; 
= (2 — 4) ((5 — A)(—4 — A) + 18) 
= (2 — A) (A? — A —2) 
= —(A — 2)*(A + 1). 

So the eigenvalues are 2 (a repeated eigenvalue) and —1. 


To find corresponding eigenvectors we solve 


—6 0 —4—),]|\w 


putting A equal to each eigenvalue in turn. 
Case 4 = —1 gives 

6u + 3w=0 

3u + 30+ 3w=0 

—6u — 3w = 0. 


Putting u = k, the first and third equations give w = —2k. The second equation 
then gives v = k. So an eigenvector corresponding to the eigenvalue —1 is 
(1 1 —2)!. 


Case 1 = 2 gives 
3u + 3w =0 
3u + 3w =0 
—6u — 6w = 0. 


So, putting w = k we have u = —k. Now v can take any value / say. Thus, any 
vector of the form [—k l k}" is an eigenvector. In particular 

ae -4 Oj’ and[—-1 0 _ 1)? are linearly independent eigenvectors 
corresponding to the eigenvalue 2. 


In this case, even though there are only two distinct eigenvalues, there are 
sufficient linearly independent eigenvectors (three) for the conditions of Procedure 
1.2 to be satisfied. Thus the general solution of the given system is 


1 0 —1 
ee Ci) tie 407 oC, Ore. 
—2 0 1 


In terms of x;, x2. and x3, the solution is 
x = Cre — Cyc" 
%a = C,e ' + C,e* 
ao —2C,e * +e”. 


Although the above example requires considerable calculation, there is little that is 
new in it. Most of it is a use of matrix methods you have met in earlier units. 


An example with complex eigenvectors 


In attempting to apply Theorem 1, it is quite possible that we may find that some 
or all of the eigenvalues and eigenvectors of the matrix B are complex. 
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Example 6 

Find the general solution of 
Xi — X32 (10) 
xX = —X 


Solution 


We have x = Bx where 


ao 


To find the eigenvalues of B, look at the characteristic equation: 
—) 1 
se |e 
ie. A7+1=0. 


So the eigenvalues are i and —i. 


To find corresponding eigenvectors we solve 


Tt all:]-° 


putting A equal to each eigenvalue in turn. 
Case A = i gives 

—lu+v=0 

—u—iv=0. 


The first equation is i times the second and so both give v = iu. Thus an eigenvector 


corresponding to the eigenvalue i is [1 iy. 
Case 4 = —i gives 
iu+ov=0 
—u+iv=0. 
The second equation is i times the first and so both give u = iv. Thus an eigenvector 
corresponding to the eigenvalue —i is [1 —i]’. 
So by Theorem 1, the general solution is given by linear combinations of x; and x, 
where 
2s ett cost + isint cost] f[sint 
Xj = : et —s eae — ; ; = ; = l . 
i ie —sint + icost —sint cos t 
and 
ie see "aol cost — isint cost] f[sint 
X2 —— ; e ut = ones = ; : == : —l . 
—j —ie™" —sint —icost —sint cos t 
Since both x, and x, are linear combinations of [cost —sint]’ =4(x, + x.) 


; 1 
and [sint cost]? = 5 1 — Xz), the general solution can equally well be written 


cos t sin t |. 
= C, : + C, . 
—sint cos t 


For our purposes this is a much better way of writing the general solution because 
we are only interested in real solutions. 


The approach in the above example works because x, and x, form a complex 
conjugate pair with real part [cost —sint]’ and imaginary part 
[sin t cost]’. It turns out that systems for which B has complex eigenvalues 


Step 1 


Step 2 


MST204 22.1 13 


can always be treated in this way provided B is real (i.e. all the elements of B are 
real). The justification for this is provided by the following theorem. 


Theorem 2 


If the conditions of Theorem 1 are satisfied, and in addition B is 
real, then the terms in the general solution 


x= C,a,e°% + Come +...+Cae™ 
a is the vector obtained by 


taking the complex conjugate 
of each element of a. 


are either real or occur in conjugate pairs ae“ and ae”. 


In Example 6 we’used Procedure 1.2 to find the general solution in the form 
x = C,ae” 5 Cae" 


where a = [1 i]’ and A =i. We then used the fact that the terms in this 
solution are conjugate to write the solution in the real form 


x = C,Re(ae“’) + C,Im(ae”) 


where Re(ae*) = [cost —sint]" is the real part of ae“ and 
Im(ae“) = [sint — cost]’ is the imaginary part of ae“. 


More generally, Theorem 2 states that if the general solution in Procedure 1.2 
contains complex terms then they occur in conjugate pairs. It is therefore possible 
to obtain a real form for the general solution by using the following procedure. 


Procedure 1.2(a): Applicable when the matrix B in Procedure 1.2 is 
real but has some complex eigenvalues 


Suppose that a is a complex eigenvector corresponding to the 


eigenvalue J. Then replace the terms ae” and ae“ appearing in 


Procedure 1.2 by Re(ae*‘) and Im(ae”). 


The method of solution of x = Bx given in this section only works if B has enough 

linearly independent eigenvectors. In applications the matrix B often turns out to The eigenvectors of symmetric 
be symmetric. When this happens Procedure 1.2 is always applicable, for ee ver discussed cee 
symmetric matrices always have a full set of linearly independent eigenvectors. at inmend ol Unit 27, Section 1. 


Exercise 6 
Find the general solution of each of the systems below. 
(1) Xi =e 6 a + 4x, 

X2= —-X1 
(ii) X,; = 5x, — 6x2 — 6x3 

X> = 4x, + 2x3 

X3 = 3X —_ 6x> = 4x3 

(Hint: 2° — 547 + 84 —4 = (A — 1) (A — 2)’). 
—2x, + 2X» 


(ili) X; 
X2 


(Solution on p. 32] 


Summary of Section 1 


We can use matrix notation to represent suitable systems of differential equations. 
A linear first-order system is in normal form if it is written 


x = Bx + k. 


We can use Procedure 1.2 to solve a homogeneous, constant-coefficient system in 
the normal form, x = Bx, so long as B has enough linearly independent 
eigenvectors. We saw also in Procedure 1.2(a) how to write the solution in an 
explicitly real form, when some of the eigenvalues and eigenvectors of B are 
complex. 
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2 Further methods for linear first-order 
systems 


2.1 Constant-coefficient inhomogeneous systems 

Suppose we wish to solve an inhomogeneous system x = Bx + h(t), such as 
Xi as 2X; + 3x, + l 

f 
X27 = 2X1 + X2 + At. 


So long as B has ‘enough’ eigenvectors (as in Theorem 1 of Section 1), it is 
possible to transform such a system into a much simpler system. We shall see 
shortly that by setting x; = 3y, + y2 and x, = 2y, — y2, System (1) becomes 


y=4yt+t 
yz = —y2 — 21. 


Each equation in (2) can be solved by the integrating factor method of Unit 2, 
since each equation now involves only one unknown variable. We can then obtain 
the solution of (1) from that for (2). 


(2) 


Before describing the general procedure, let us see an example of how it works. I 
shall make a note of the steps involved in the margin. 


Example 1 
Solve System (1) above. 


Method 
The system may be written 


x = Bx + h(t) (3) 
2 3 . 
where B = 5 ' and h(t) = [t 4t]". 


The first step is to find the eigenvalues and eigenvectors of the matrix B. In fact we 
have already done this in the solution to Exercise 5 of Section 1 where we found 


that: 
[3 2] is an eigenvector corresponding to the eigenvalue 4; 
[1 —1]’ is an eigenvector corresponding to the eigenvalue —1. 
The next step is to form a matrix P, whose columns are the eigenvectors. So here, Are there enough 
| eigenvectors? If so, form P so 
P that P~ ‘BP is diagonal. 
P= 5 it 


We then know (by Theorem 5 of Unit 21, Section 3) that P is non-singular and 
that 


0 
P~'BP = ! 


That is, P~'BP is a diagonal matrix with diagonal entries equal to the eigenvalues. 
We now substitute x = Py in (3), where y is a new vector of variables. This gives Put x = Py 


Py = BPy + h(t). 


(Here we have used the fact that since all the entries in P are constant, 
differentiating Py with respect to t gives Py.) 


Since P is non-singular we can multiply by P~' to obtain 
y=P 'BPy + P ‘h(?). Calculate P~'. 


Because P is a 2 x 2 matrix we can easily calculate its inverse: 


ees ee 
pote | 
5 : 
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Thus, in terms of the y variables, System (1) becomes 


ee a eae 


That is, 
eee eo 
yo = —y2 — 2. 


We now have a system that we can solve, because each equation can be solved by 
the integrating factor method of Unit 2. For example, the first equation can be 
written 


Yi —4y1 = 0. 
Multiplying both sides by the integrating factor e-!*“' = e~* gives 
e- “(1 — 4y1) = te 
that is 
d 4t + 
=e 9 = le". 
dt (yie “) = te 


Integrating both sides we obtain 
SPY ies = re-*ae a C, 


=(-it-zy)e “4+ C, (integrating by parts), 
thus 
jy, = Ce" —ft— 
The second equation may be solved similarly to give (omitting the details) 
yz = Cyne ' — 2t + 2. 
Now the solution of the original system is x = Py. That is, 
= . } 2 pst —4t— a 
X> 2 —1]| Coe"' —2t+2] 
Thus 
Xx, = 3C,e“ + Coe '-—Yt+# 
X_ = 2C,e" — Cyne ' 4+ 3t — 4. 
So we have solved the original system (1). 


(The result of such a calculation should of course be checked, though I will not 
give details here.) 


We could have used a method similar to the above for the homogeneous systems 
in Subsection 1.2, but there the above method would have been unnecessarily 
complex. Here we need to make the transformation x = Py to get the system into 
a ‘diagonal’ form, that can be solved equation-by-equation using the integrating 
factor method. 


Solve the resulting equations 


48 


y=P ‘BPy+P'h 
individually, by the integrating 
factor method. 


Now find x, from x = Py. 


Check 
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The procedure for inhomogeneous systems is set out below. 


Procedure 2.1: To solve a linear, constant-coefficient, 
inhomogeneous, first-order system 


For this procedure to work we need to be able to express the 
system in the normal form 


x = Bx + h(t) 


where B is an n x n matrix with n linearly independent 
eigenvectors. 


1. (i) Find n linearly independent eigenvectors of B; 
(11) Form a matrix P whose columns are these eigenvectors; 
(iii) Find P™?. 
Put 
x = Py 
in (3), to obtain 
y =P 'BPy + P ‘h(t), (4) 
where P~'BP is diagonal. 


Find y by solving the equations in (4) individually, using the 
integrating factor method. 


Transform the solution for y back to a solution for x, using 
x = Py. 


Exercise 1 

Solve the system 
X, = 2x, + 3x, +e” 
X_ = 2x1 + X2 + 4e”. 


[Solution on p. 34] 


Exercise 2 
Solve the system 


x; = —2x, i 10x, + 48 
X, = 10x, — 2x, + 9le’. 
[Solution on p. 34] 


2.2 Numerical methods 


We saw in Units 2 and 6 that differential equations cannot always be solved by 
analytical methods. In such a situation we may resort to obtaining numerical 
approximations to the solution. Techniques for doing this were discussed 
extensively in Unit 19. 


We are now faced with a similar situation for systems of differential equations. The 
analytical techniques described in this unit are only applicable to certain special 
types of systems and so the need to have numerical techniques for systems of 
differential equations is as pressing as it was for single differential equations. 


Although, for reasons of space, I cannot give a comprehensive treatment of the 

topic, it is useful to indicate here how we can go about constructing numerical 

schemes to generate approximate solutions to first-order systems of differential Euler’s method was discussed 
equations. We will look at a numerical method based on Euler’s method for a in Units 2 and 19. 
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single differential equation. Recall that in Euler’s method for a single differential 
equation we approximate the solution of the initial condition problem 


x 
ae = nd ¥ = X_ when t=, 


by starting with Xo = Xo and using the recurrence relation 
X41 =X, + hm(X,,t,) | 

to generate approximations X, to the true solution x(t,) at t, = to + rh. 

We can easily extend Euler’s method to a system of two equations of the form 
Xy = M4 (X1,X2, Ct) 
X2 = Mz(X1,X2, ¢). 


We do this by using Euler’s method to generate approximate solutions to each of 
these equations. Looking at the first, Euler’s method gives an approximation X,, 
to x,(t,) where 


P, ere - Aes + hm, (X15, As 5s t,). 


Here X,, is the approximation to x2(t,) where (applying Euler’s method to the 
second equation) 


oo —— X 2» + hm, LX e5. X29 t,). 
This may seem a little complicated but in practice it is quite straightforward. Let 


us look at an example. 


Example 2 


Write down the recurrence relations corresponding to an application of Euler’s 
method to the system 


X, = tx; -—tx, 
X, = 3x, — cost, 
Where x, = 1.x. = -—3 ai=@ 
Solution 
Applying Euler’s method to each equation in turn gives 
Xieet =X1, + h(tX1, — X2,) 
Xort1 = X2, + h(3X,, — cost,). 
We have also that X; 5 = 1, X29 = —3, and that t, = rh. 


(5) 


The recurrence relations (5) obtained in the Example can be used to generate 
approximations to the solution of the given system, step-by-step. Starting from the 
known values of X; 9 and X20, we use (5) to calculate X, , and X2,. We then 
use (5) again to calculate X,, and X,; and so on. The results obtained depend 
on the step-length h chosen, and will usually be more accurate the smaller h is. 


Example 3 


Use the recurrence relations (5) to calculate X;,. and X2. with h = 0.1. What 
values of the solutions do these approximate? 


Solution 
Putting r = 0 in (5), we obtain 
X11 = X19 +01 (toX1,.9 — 18X20) 
X21 = X29 + 0.1 (3X 1,0 — COS to). 
Substituting X; 9 = 1, X29 = —3, and to = 0, we obtain 
Xi,=1 
X2.,= —3+0.1(3 —1) = —-2.8. 
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Putting r = 1 1n (5) gives 
Xia = X11 +01 O1X,, 001%, J 
=1+40.1 0.1 + 0.028) 
= 1.0128. 
X22 = —2.8+40.1 (3 — cos0.1) 
= —2,.599S. 


Since we are working with step length h = 0.1, X 1,2 1s the calculated 
approximation to x,(0.2), and X, > the approximation to x (0.2). 


Calculations such as this are time-consuming and would normally be performed 
on a computer, so I shall do no more by hand here. 


It is a straightforward matter to extend the above method to systems containing 
more than two equations. We can consider systems of equations which can be 
written in the form 


X4 se MANTA hes t) 


X2 sd M2(X1,X2,...,Xny t) 


Xn = M,(X1,X2, veo Xny t). 


If we introduce the vector notation x = [x, X> coe x,]' these equations 
can be written more concisely in the form 


x = m(x, f). 


The procedure for applying Euler’s method to such systems is given below. 


Procedure 2.2: Euler’s method for systems of equations 


This procedure applies to systems of the form 


xX = m(x, f) 


with initial condition x = xo at t = fo. 


1. Replace each equation in the system by a recurrence relation, 
using Euler’s method. Thus the equation 


X; = m;(x, t) 


is replaced by 
| Fare ae Xir - hm;(X,, iJ. (i = ee (8) 


Here X= [X, . X33 ae Xn7]' is the 
approximation to x(t,) where t, = to + rh. 


EO SA%, = x, 


(ii) Choose a step length h. 


(ii1) Use the recurrence relations (8) to calculate X,, then use 
(8) again to calculate X,, and so on. 


Exercise 3 


(i) Write down, in a form suitable for generating solutions, the recurrence relations 
corresponding to an application of Euler’s method to 


xs = 3tx +4 
X2 = tx; —x,—-—4?, 


where x, = 5,x, =2att=0. 


MST204 22.3 19 


(ii) Use the recurrence relations you obtain, with step-length 0.1, to calculate 
approximations to x,(0.2) and x,(0.2). 


[Solution on p. 34] 


We have seen in previous units that care must be taken in the use of numerical 

methods. There is no time here for a discussion of the important topic of stability, For a discussion of stability 
which is essential to a proper understanding of numerical schemes such as those | see Unit 19. 

constructed above. Approaches other than those based on Euler’s method are 

possible; in fact we could develop methods for first-order systems based on any of 

the methods for a single first-order equation described in Unit 19, but, again, I will 

not pursue this topic here. 


Summary of Section 2 


To solve an inhomogeneous constant-coefficient system of differential equations 
x = Bx + h(t) where B is an n x n matrix with n linearly independent eigenvectors, 
we use Procedure 2.1. 


A method of constructing recurrence relations to generate approximate solutions 
to systems of differential equations was described. This is based on Euler’s method 
and consists of applying the approximation specified by Euler’s method to each 
equation in the system in turn. 
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3 Second-order homogeneous systems of the 
form Xx = Bx: simple harmonic motion 


We now turn our attention to systems of differential equations which involve 
second derivatives. We shall be concerned solely with linear, constant-coefficient, 
second-order systems. Such systems arise frequently in modelling, especially in 
mechanics. A particularly important class of second-order systems are those that 
can be written in the form X = Bx, and we study these in this section. 


Subsection 3.1 summarizes material covered in the television programme ‘Two-way 
motion’. The programme looks at the solution of a particular system of the form 

x = Bx and illustrates the sort of real situation in which such a system of 
differential equations may arise. In Subsection 3.2 we look at a procedure for 
solving systems of the form xX = Bx. 


3.1 ‘Two-way motion’ (Television Subsection) 


Pre-programme notes 
In the programme we study the solution of the pair of differential equations 


X = —3x + $y 

See (1) 

A meee mage =e 
This pair of differential equations provides an approximate description of the 
motion of a ball-bearing in a suitably shaped bowl (where x and y are the 
horizontal position coordinates of the ball-bearing). The programme consists of 
three parts: 
1. The actual motion of the ball-bearing in such a bowl; 
2. A computer animation of the precise motion implied by Equations (1); 
3. The mathematical solution of Equations (1). Figure 1. We consider the 

: motion of a ball-bearing in a 
(The computer animation enables us to see how accurately Equations (1) describe bowl. The two horizontal _ 
the actual motion of the ball-bearing. For various reasons the actual motion is components of acceleration, x 
: : and jy, are specified 

only rather crudely approximated by Equations (1). For example the equations (approximately) by the 


take no account of various factors tending to damp the motion.) Equations (1). 
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To solve the equations mathematically it is helpful to use vector notation. If we 
write r for the vector of variables [x y]’ and A for the matrix 


5 4 
4 5 
5.9 
then we see that the pair of differential equations (1) can be written 
f= Ag. 
This equation is similar in form to the single differential equation 


j= ky. 


For k < 0, this is the equation of simple harmonic motion. The similarity of the 
equations Ff = Ar and jj = ky is not simply a matter of mathematical formalism. 
For suitable matrices A, the vector equation Ff = Ar does describe a generalization 
of simple harmonic motion to motion in more than one dimension. The modelling 
of mechanical systems displaying such motion is discussed in Unit 24 on normal 
modes. Our concern here is with the mathematical solution of the equation Ff = Ar. 
However, in the programme we take an informal look at an example of two- 
dimensional motion—the ball in the bowl. This example is provided to help you 
to visualize the solution of the equation Ff = Ar which is obtained in the 
programme. 


The solution of the equation involves knowledge of the eigenvalues and 
eigenvectors of the matrix A. I will ask you to find these before viewing the 
programme. 


Exercise 1 
Calculate the eigenvalues and eigenvectors of 


[Solution on p. 35 | 


Now view the television programme ‘Two-way motion’. 


Post-programme notes 
To solve a system of the form 


¥ = Br (2) 
we may look for solutions of the form 
r= ee" (3) 
where a is a constant vector, and a # 0. With r given by Equation (3), we have 
f=pae" and = *F= p’ae 
Hence (3) is a solution of f = Br if, and only if, 
p’ae" = Bae: 
i.e. u?a = Ba. 


With a # 0, this condition asserts that yu? is an eigenvalue of the matrix B, and a 
is a corresponding eigenvector. 


If A is an eigenvalue of B and a is a corresponding eigenvector, we have two 
solution of f = Br of the form r = ae“. These are 


r = ae*'and r = ae~V44. 


In the programme, we concentrated on the pair of Equations (1) on page 19, 
which are equivalent tof = Ar withr= [x _ y]’ and 


ah 4 
a-| | (4 


You will find this notation 
discussed further in Subsection 
| SS 


Simple harmonic motion was 
discussed in Unit 7. 


The ball actually moves in 
three directions but its motion 
is completely determined by 
the behaviour of its two 
horizontal coordinates x and y. 


TV 22 


We are not using exactly the 
same notation as in the 
programme. 
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The eigenvalues and eigenvectors of A are given in the solution of Exercise 1, 
where we found that: 


[1 —1]’ is an eigenvector corresponding to the eigenvalue —1; 
[1 1]’ is an eigenvector corresponding to the eigenvalue —3. 


So, the solutions of F = Ar given by the eigenvalue —1 are 


a = er 
| if and| _| | 


More generally, a linear combination of these solutions is also a solution. That is 


1 . 
i E fee + d,e™"*), 


where c, and d, are constants, is a solution. Use of Euler’s formula enables us to 
rewrite this as 


= |i |e cost + D, sin t) (5) 


where C, and D, are arbitrary real constants. Similarly, corresponding to the 
eigenvalue —%, we have a solution 


r= | fese™ * d,e— Stes! 
1 


or 
1 1 ae | 
i= \ (C,cos3t + D,sin3f). (6) 


Although (5) and (6) contain arbitrary constants, neither is the general solution of 
r = Ar. The general solution is a linear combination of both (5) and (6). That is 


1 1 
= aCe cost + D, sint) + he cos 3t + D, sin 3t). (7) 


The programme showed computer animations of the motion of a point in a plane, 
where the position coordinates (x, y) of the point satisfy the original pair of 
differential equations (1) (or, equivalently, Ff = Ar with A given by (4) and 

r=[x  y)’). | 

Each of the Solutions (5) and (6) above corresponds to a particularly interesting 
motion of such a point. Looking at (5) we see that it gives the position 

y= [x y]’ asa scalar multiple (C, cost + D, sint) of a fixed vector [1 —1]’. 
So, in the motion described by (5), the point moves along a straight line defined 
by the vector [1 —1]? (see Figure 2). The factor C, cost + D, sint varies 
sinusoidally, so that the point oscillates in simple harmonic motion along this line. 


Similarly Equation (6) implies simple harmonic motion along the line defined by 
the vector [1 1]’. Notice also that Equations (5) and (6) show that the 
angular frequencies of these two special motions are in the ratio 1:3. 


The motions described by Equations (5) and (6) are only achieved under suitable 
initial conditions. The general motion of the point, under arbitrary initial 
conditions, is given by the general solution (7) of the equation Ff = Ar. But any 
solution of the form (7) is just the sum of a solution of the form (5) and one of the 
form (6). Thus any motion of the point can be seen as the sum of two special 
motions, one along each of the lines defined by [1 1}’ and [1 —1]’. The 
programme included a computer animation illustrating this fact. Figure 3 shows 
examples of the paths (orbits) that such a point may follow. The general motion 
described by the Equations (1), or their solution (7), is of the point tracing and 
retracing orbits of the form shown in Figure 3. 


Now consider the real bowl. By viewing the bowl from above we observe the 
behaviour of the horizontal x, y co-ordinates of the ball bearing. By comparing 


Euler’s formula was discussed 
in Unit 5. It states that 


e”® — cos@ + isin 0. 


Figure 2. If the ball-bearing is 
released from either of the 
positions A or B, then the 
Equations (1) imply simple 
harmonic oscillations along 
the straight lines /, or /, 
respectively. Such straight line 
motions were observed in the 
real bowl. 


(a) Motion of the point 
starting from rest at A. 


(b) Motion of the point 
starting from B in the 
direction shown. 


Figure 3. Examples of orbits 
determined by the Equations (1). 
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this view of the bowl with the computer animations we can judge how accurately 
Equations (1) describe the motion of the ball bearing. In the programme we found 
that the ball bearing could be made to describe oscillation along the lines defined 
by [1 1]’ and [1 —1]’. However, the oscillations were damped rather than 
pure simple harmonic motion. Furthermore, the general orbits of the ball bearing 
in the real bowl were not of the form shown in Figure 3, for as well as damping, \\ 
there were other inaccuracies in the real bowl. One was that the frequencies of the d 
two special motions were not exactly in the ratio 1:3. We concluded the 
programme by trying to adjust the mathematical solution so as to make it give a 
more accurate description of the behaviour of the real ball-bearing. We did this by 
looking at computer animations obtained by adding simple harmonic motions in 
the same two directions [1 . 1]? and[1 —1]’ but with frequencies in a ratio 
which differed slightly from 1:3 (see Figure 4). 
Figure 4. By adjusting the 
Exercise 2 ratio of the animation’s 


| ss a tT: . ; frequencies in the special 
An object moves so that its position, r = [x y]’, in a plane satisfies the equations directions, we cin obtain a 


X = —4Px + Sy better correspondence between 
the animation and the motion 
¥ = 9x — Py. in the real bowl. 


Give the directions (if any) in which such an object can describe simple harmonic motion in 
a straight line, and the frequency of any such motion. 


[Solution on p. 35] 


3.2 Solution of second-order systems of the form X = Bx 
In the previous subsection we saw that 

x= ae” (4) 
is a solution of 

x = Bx (3) 


so long as ? is an eigenvalue of B and a is a corresponding eigenvector. For each 
non-zero eigenvalue 4 of B, we have two solutions of (3): aev*' and ae~¥*". Aside 
from the fact that there are more terms in the solution, solving a second-order 
system of the form (3) is no harder than solving a first-order system. The general 
solution of (3) can usually be built up from solutions of the form (4). But as in 
Subsection 1.2 this method does not invariably work—it requires that the matrix 
B has enough eigenvectors. The following theorem and procedure for second-order 
systems of the form X = Bx are very similar to those in Subsection 1.2 for first- 
order systems. 


Theorem 1 


Let B be an n x n matrix. Suppose that B has n linearly independent 
eigenvectors a;,4>,...,a,, corresponding to the eigenvalues 
Ay,A2,..+, Any Tespectively. 


If all the eigenvalues are non-zero the general solution of the 
system of differential equations 


X = Bx 


n pattie 
Sey alee ae. 
r=1 


where c, and d, are arbitrary constants (1 <r < n). 


If any eigenvalue 4, is zero the general solution is obtained by 
replacing the rth term in (5) by (c, + d,t). 


In the example studied in the previous subsection the eigenvalues were negative 
and we had to express the general solution in terms of real functions. In general, 
the real form of the terms c,ev*"' + d,e~¥*"' in (5) depend on the nature of the 
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eigenvalues. There are three cases to consider: (1) A, real and positive, (2) A, real 
and negative and (3) 4, complex. (The case 2, = 0 is considered separately in 
Theorem 1.) 
Case 1: A, real and positive 
In this case the square roots of 4, are real and the term in (5), 

aioe de, 
is already explicitly real. 


Case 2: i, real and negative 
Now the square roots of 4, can be written as 
i./—A, and —i,/—A,. 


Since both e'¥~*' and e~'¥~*' are combinations of cos J —A, t and sin./—A, t, 
the term in (5) can in this case be written as 


a,(C,cos./—A, t + D, sin ./ —A, t). 
where C, and D, are combinations of c, and d, and are again arbitrary constants. 


Case 3: 1, complex 


In this case we can replace the complex exponentials in (5) in a way similar to that 
used in Procedure 1.2(a) for homogeneous first-order systems. 


We know from Theorem 1 that if 2 is a complex eigenvalue of B with 
corresponding eigenvector a then the functions 


aev*' and ae~v?' (6) 
will appear in the general solution (5). It turns out that if the matrix B is real then 
4 will also be an eigenvalue of B, and 4 will be a corresponding eigenvector. This 
implies that the functions 

aev*' and ae~v*' (7) 
will also appear in the general solution. Now the functions (7) are complex 
conjugates of the functions (6) and so we can replace all four functions by the real 
functions 

Re(aev*"), — Im(aev*"), — Re(ae~¥*"), = and_~—sIm(ae~¥“%), 

To summarize, we have the following procedure: 


Procedure 3.2: 
To solve a linear, second-order, system of the form X = Bx 


For this procedure to work B must be an n x n matrix with n 
linearly independent eigenvectors. 


1. Find the eigenvalues 4,,/2,...,4, of B and a corresponding 
set of linearly independent eigenvectors a,,a>,...,a,,. 


Write down the general solution in real form. There are four 
cases to consider: 


(i) If all the eigenvalues are real and positive the general 
solution is 


n 


x= a(C,evert + D,e~¥?*r?). (8) 


r=1 


If A, = 0 replace the rth term in (8) by 
a.(C, + D,t). 
If A, is real and negative replace the rth term in (8) by 


a,(C, cos ./ —A, t + D,sin./ —A, t). 


If an eigenvalue / is complex replace the functions aev**, 
ae~¥*' gev"' and ae~Y*' in (8) by the functions 


Re(aev* a: Im(aev* 2 Re(ae~v* ‘Yand Im (ae~v* y. 
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The eigenvalues of the matrix B do not necessarily all fall into the same case, of 
course. The example below illustrates this. 


Example 1 
Find the general solution of the system 


X¥, = X1+3x,- X3 

X. = 3x, + X2+ X3 

x3 = 4x3. 
Solution 


We have a system of the form X = Bx with 


1 bE | 
B= |3 1 re. 
0 0 4 


The first step is to find the eigenvalues and eigenvectors of B. 


The eigenvalues can be found by solving the characteristic equation 
det(B — AI) = 0. We have 


i Se 
dt 1H =|. SS 


=. 
1-—A 


= (4— 4(l — af -— 3) 
= —(4—A)(4—A)(2+A4). 


—) 
= 4-0); 


So the eigenvalues are —2 and 4 (repeated). 


The eigenvectors corresponding to the eigenvalue —2 satisfy 


e 


3 3 —Illju 0 
3 3. Tt fe j=] 0}. 
0 0 6} \w 0 
Thus w = 0 and u = —», and so an eigenvector is [1 —1 0]’. The 


eigenvectors corresponding to the eigenvalue 4 satisfy 


—3 3 —Illlu 0 
3 -—3 (ea aie ew 6 is 
0 0 0} |w 0 


These equations reduce to 
—3u+3v—w=0. 


Since we effectively have only one equation in three unknowns, we can find two 
linearly independent eigenvectors corresponding to the eigenvalue 4: for example 
(1 1 OF aa © SI a 


By Procedure 3.2 the general solution is 


1 0 1 
“=| 1 KC e* + Die *) +11 (C,e" + Doe 1 pcos </21t + Disin 2/21) 
0 3 0 
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or alternatively 
x, = Ce" + Die * + C3cos,/2¢ + D3sin,/2t 
x2 =(C, + Cre" + (D, £ De" *= Cytos/ 21 —~ D3sin,/2t 
Xs = 3(C,e" + De *), 


Exercise 3 
Find the general solution of each of the systems below. 
(i) X,;=x, 
X> =X 
(ii) x — — 10x, em 6x> 
2X> = —12x, —_ 20x>. 
[Solution on p. 35] 


Summary of Section 3 


To solve a second-order, constant-coefficient, system of the form X = Bx, for 
suitable B we may proceed in a similar way to that described in Section 1 for a 


first-order system. The details are given in Procedure 3.2. 
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Of particular significance in mechanics is the case where the eigenvectors of B are 


real and negative. In this case the model ¥ = Bx represents a generalization of 

simple harmonic motion (Unit 7). In the television programme we looked at an 
example of this: the two-dimensional motion of an object whose position vector 
r= [x y]’ varied with time according to the system of differential equations 


r = Ar where the matrix A had negative eigenvalues. 


4 Forced oscillations 


4.1 The structure of the solution 


In the unit ‘Differential equations II’ we introduced certain results about the 


solution of the forced oscillation equation 


X + 20x + w?x = b, cos vt + bz sin vt. 


These results were used extensively in the mechanics unit on ‘Damped and forced 


vibrations’. 


In this section, we discuss the extension of these results to systems of differential 


equations. I will first remind you of the results for a single equation, by looking at 


an example. 


Example 1 
We consider the solution of the equation 


X + 2x + 5x = 20083t. 


To solve this equation we first calculate the complementary function; that is, the 


general solution of the associated homogeneous equation 
eo oy + Sx = t. 
I will not give details of this calculation. The complementary function is 


xX, =e '(Acos2t + Bsin 2t). 


(1) 


Next we calculate a particular solution. We can do this by the method of phasors, 


in which we try x = Re(ze*") and calculate the complex number z. Since the 


phasor of 2 cos 3t is 2, we get here 


((3i)? + 2(3i) + 5)z = 2, 


and so 
2 


6i — 4 


| 


Note the 2 on the left-hand 
side of part (ii). 


Unit 6, Section 4 


Unit 8 
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Thus a particular solution is 
Xp = —73C0s 3t + #ssIin 3¢. 


The general solution of (1) is found by adding the particular solution and the 
complementary function, and is 

X=X,+Xy, 

=e '(Acos2t + Bsin 2t) — ~;cos 3t + sin 3¢. 

Now let us consider the behaviour of this solution for large values of t. The 
complementary function tends to zero for large t, because of the term e~‘'. Thus for 
large t any solution of (1) is approximately equal to the particular solution 

Xp = —73008 3t + Fy sin 3¢. 


This is known as the steady-state solution of (1). We refer to the complementary 
function as a’ transient, since its effects die out as t becomes large. We saw in Unit 8 
that in many mechanical problems on forced vibrations, a knowledge of the 
steady-state solution is all that is needed, so long as we know that the 
complementary function is indeed a transient (i.e. it tends to zero for large t). 


In general, we have the following important information about the solution of the 
single equation 


A,X + a,X + a3x = b, sinwt + b, cos at. (2) 


(1) The general solution of (2) is the sum of the complementary function and 
any one particular solution. 


(11) So long:as each of a;, az and a; is positive, the complementary function is a 
transient (i.e. it tends to zero for large f). 


(111) There exists one purely sinusoidal particular solution, to which any solution 
converges for large tf. 


Now the motion of a system of vibrating particles under an external force may be 
modelled by a system of differential equations, of the form 


A,X + A,x + A3x = b,; cosa@t + b, sin at. (3) 


So analysis for this system parallel to that above for the single equation (2) would 
be very informative. To achieve this we need three things. 


(1) To know that the general solution of (3) is made up of ‘particular solution 
plus complementary function’. 


(11) To have a condition on the matrices A,, A, and A; which ensures that the 
complementary function is a transient. 


(111) To know that there is a purely sinusoidal solution of (3), and how to 
calculate it. 


Each of these can be dealt with, and we will do so now. Points (i) and (ii) are 
resolved by Theorems 1 and.2 below. These are stated without proof, although a 
proof of Theorem 1 would be quite straightforward. The proof of Theorem 2 
depends on methods in linear algebra beyond the scope of this course. 


Theorem 1 


Suppose that x, is a particular solution of the system 


A,x i A>x + A3X = h(t). 
Then the general solution of (4) is 
X=X),+ X, 


where x, (the complementary function) is the general solution of 
the associated homogeneous system 


A,X Si A>x + A3X = 0. 


We shall study such systems 
in Unit 24. 
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Theorem 2 


Suppose that each of the matrices A,, A, and A; has constant 
real entries, is symmetric, and has the property that each of its 
eigenvalues is positive. Then any solution of the system 


A,x + A>x + A3x = 0 


has the property that x tends to 0 as t becomes large. (That is, 
every component of the vector x becomes small as t becomes 
large.) 


Exercise 1 


Which of the systems below satisfies the criteria given in Theorem 2 for the solution to be a 
transient? 


(i) a + 3x, + 2X, + 2x, —xX,=0 
3X2 + 2X, + 3X, — x, + 2x, =0 


(ii) X, + 3%, + 2X2 + 2x, =0 


Gti) Xy +X, = 2x, + x5 
Xz +X2 =X, + 2x. 
[Solution on p. 36] 
Exercise 2 
A particular solution of 
X,=X2,+1 
X2=x,+1 


is not difficult to spot. Do so, and hence find the general solution of the system. (You may 
reuse the solution of Exercise 3(i) of Section 3 here.) 


[Soluton on p. 36] 


4.2 Finding a sinusoidal particular solution 


In (ii) above, we wished to know a purely sinusoidal solution of a system of 
equations of the form (3). In this subsection we shall see how such a solution can 
be calculated. The method uses phasors and is an extension of the procedure in 
Unit 6 for finding a sinusoidal particular solution of a single second-order 
differential equation. We look for a solution of the form x = Re(ze'“‘) where z is a 
constant vector with complex entries. Let us see an example. 


Example 2 
Find a sinusoidal particular solution of the system 
X; +X, =sin2t 
X2 +X, + 2x, =cos2t. 
Solution 
Step 1: First we write the right-hand sides in phasor form: 
X; +X. = Re(—ie?*) 
X. +X, + 2x. = Re(e?"). ” 
Step 2: Next we look for a solution of the form 
x = Re(ze?"*), 
Then if x = [x, x2]" and z= [z, z>]', we have 
x, = Re(z,e7"), 
SO ¥ = ReQiz,e*"), 
and = xX, = Re(—4z,e7"); 


and similarly for x). 


A symmetric matrix A is one 
for which A? = A (see Unit 20, 
Subsection 4.4). 


Phasors were defined in Unit 5 
and used in Units 6 and 8. 
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Thus to satisfy (5) we need 


—47, + 2iz, = —i (6) 

—4z, + 2iz; + 2z, = 1. (7) 
Equation (7) may be rewritten as 

2iz,; — 2z, = 1. (8) 


Step 3: Thus we have a pair of simultaneous equations, which may be solved by 
Gaussian elimination, as follows. Adding 2i/4 times (6) to (8) gives 


—3z, =} 
i.e. Zz, = —4. 
| Substituting in (6) gives 
—47,-i=-i 
be. Zoe MD, 


Step 4: Now z, and z, are the phasors of the required particular solution. The 
solution is therefore 


X14 = 0 
X2 = —4$c0s 2t. 


The method used in the example is set out in the following procedure. 


Procedure 4.2 
To find a sinusoidal particular solution of 
A,X + A2x + A3x = b, cos wt + bz sin wt. 
Write the system in the form 
A,X + A.x + A3x = Re(de’”). 


Look for a solution of the form 
x = Re(ze'), 


where z = [2, 2, ee oe ee ez 38-8 Complex 
number. This leads to a system of simultaneous linear 
equations for the z,’s. 


Solve this system of linear equations, and so find z. 
The required solution is then 


x = Re(ze'“). 


The only significant difference between this procedure and that in Unit 6 for a 
single equation is that here, in Step 3, we must solve a set of simultaneous linear 
equations, rather than just a single equation. 


Let us see what happens if we apply Procedure 4.2 in general. Suppose we have 
written our system in the form 


A,X + Ajx + A3x = Re(de’”) 


where d is some constant, complex, vector. We look for a solution of the form 
x = Re(ze''). Then x = Re(iwze'”’) and X = Re(—w7ze'"). So we have a solution 
so long as 


A, (—w?ze'®') + Az (iwmze’) + A3(ze) = de, 
a (—w’A, s&s iwA, + A3)Z = d. 


This represents a set of n simultaneous equations (with complex coefficients) in n 
unknowns (compare with Equations (6) and (7) in Example 2). These equations 
will usually have a unique solution for z, and so we will obtain a sinusoidal 
particular solution x = Re(ze‘“”). 
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There could be problems with this method if the matrix —w?A, + iwA, + Az Is 
singular. (This compares with the problem of finding a particular solution of, for 
example, the single equation 


xX + 4x = sin 21, 


which has no solution of the form x = Re(ze?"").) So, once again, we have a 
procedure that will usually work, but which could break down. I will not get into 
a discussion of how to deal with possible ‘problem cases’ here. 


The question below involves the use of Procedure 4.2. You may wish to treat it as 
an Exercise, or as another worked Example, before going on to try Exercise 3 for 
yourself. 
Question 
Find a sinusoidal solution of 

X, + 8X2+x, = 8sin3t — 8cos3t 


Xo+ X = —6sin 3t — 3cos3t. 


Solution 
Step 1: First we write the right-hand side in the appropriate form 


x, a 8X, - xX, = Re(— 8ie** 4 8e>'*) oe Re ((—8i ae 8)e*"") 


X2 +X, = Re ((6i — 3)e*"). 
Step 2: Next we look for a solution of the form x = Re(ze*"), where | 
x = [x, x2]? andz= [z, z2]'. This is a solution if 

(3i)*z, + 8(3i)z, +z, = —8i--8 

(3i)?z2 + 3iz, = 6i — 3. 
That is, if 


—8z, + 24iz, = —8i — 8 
3iz; — 92, = 6i — 3. 

Step 3: Multiplying the first equation by 3i/8, and adding this to the second 
equation, gives 

—18z, = 31. | 
i.e. z, = —#i. 
Substituting this into the first equation gives 

—8z,+4= —-8-8. 
Le: 73 =3+i 
Step 4: Now z, and z, are the phasors of the required particular solution. The 
particular solution is therefore 

x, = cos 3t — sin 3t 

Xx, = ¢sin 3¢. 


(As usual, this solution should be checked by substitution in the original 
equations. ) 


Exercise 3 
Find a sinusoidal particular solution of each of the systems below. 
(i) x; + 4x, + 2x, = 6cos 2t 
X2 ec ae 9x, = 2sin 2t 
(11) XX; +2X,.+ X,+ x, — 3x, =sint 
3X; + X2+2X,+2x,+ x, =cost —2sint. 
[Solution on p. 36] 
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Summary of Section 4 
We discussed results about the solution of the system 
A,X + Ax + A3x = b, coswt + b, sin wt (3) 


relevant to the motion of a system of vibrating particles. The results are similar to 
those for a single second-order equation. The general solution of (3) is the sum of 
any one particular solution and the complementary function: that is, the general 
solution of A,X + A,x + A3x = 0. We saw that if each of the matrices A,, A, and 
A; has constant real entries, is symmetric, and has all its eigenvectors positive, 
then the complementary function is a transient; that is, all the entries in the vector 
x tend to zero for large t. We saw also that (3) has a particular solution of the 
form x = d, cos wt + d, sin qt. This can be calculated by looking for a solution of 
(3) of the form x = Re(ze'”) and calculating the complex vector z from the 
resulting linear equations, as described in Procedure 4.2. 


Thus, so long as the matrices A,, A, and A; are of the appropriate form, any 
solution of (3) settles down to the same steady-state solution in the long-term. This 
is independent of the initial conditions, which only affect the transient term. 


5 End of unit exercises and problems 


Exercises 1-3 below cover the main objectives of the unit. They can be used either 
for further practice, or when revising as a check that you have covered the central 
ideas of the unit. 


Problems 1 and 2 take the ideas in the unit a little further. 


Exercise 1 
Find the general solution of each of the systems below. 


(i) X41 = 8x4 is 5x2 (v) xy = 8x; = 5X2 + sin 2t 
x2 = 10x, — 7x, X2 = 10x, — 7x2 + 2cos2t 

(it) xX, = 8x, — 5x, + 2e (vi) X, = —x, + 2x, 
Xo = 10x; — 7x, — e” X, = —x, — 3x, 

(iii) %s = 6x, =e 4x> ag AG 4x3 
X> = 10x, — 6x> (vil) xX; = 4x, + 4x, —_ X2 

(iv) x, = 8x, — 5x, 3x1 + X2 = —X3 + X2 

| ¥, = 10x, — 7x, Es sight 


[Solution on p. 36] 


Exercise 2 
Describe the behaviour of the solutions of the following system for large values of t. 
4x, + X2+ 3X, + 3x, + 4x, =cos2t 
X, + 4X2 + 3X, + 4x, + 6x, = sin 2t. 
[Solution on p. 39] 


Exercise 3 


Construct a set of recurrence relations, based on Euler’s method, suitable for generating 
approximations to the solution of the following system. 


Xi =X. 

X2=X3 

X3 = xX, + 3x2 — tx3 — cost, 
where x,(0) = 2, x,.(0)=-1, x3(0) =4. 
[Solution on p. 39] 
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Problem 1 
(i) Suppose that y = x, and that x, satisfies the system given in Exercise 3. Show that y 
satisfies the differential equation 
y+ ty —3y —y = —cost. 
(ii) Can the process in (i) be reversed? Given that y satisfies the differential equation 
V+ 4ty = sin 2t, 
with y(O) = 1, y(O) = 2, (0) = 0, can you construct a first-order system that y must 
satisfy ? | 
Use this idea to construct a set of recurrence relations to generate numerical solutions 
to this differential equation. 


(iii) The idea introduced in (ii) above can be extended to systems of higher-order equations. 
Consider, for example, the pair of second-order (non-linear) equations below, which 
model the motion of a particle in two dimensions under the influence of gravity, and of 
a resisting force whose size is proportional to the square of the velocity. 


i — A(x? + y7)"?x 


x? + y? 
i= A? + yy, 
ae - 
Here k and J are constants. 


By introducing new variables x, = x, x, = y, x3 =X, X4= y, construct a first-order 
system equivalent to the second-order system. | 


Hence construct a set of recurrence relations that would generate numerical solutions 
to this second-order system (for given values of k and 4). 


[Solution on p. 40] 


Problem 2 (This introduces methods that cover the special cases not dealt with by the 
Procedures in Sections 1 and 2.) 


(i) Show that x = ae“ + bte* is a solution of the system x = Bx so long as 

(a)b is an eigenvector of B corresponding to the eigenvalue 4, and 

(b)a satisfies (B — Al)a = b. 
(11) (a)Try to solve the system 

Xs = Xs 
X, = —4x, + 4x, 
by Procedure 1.2 of Section 1. What goes wrong? 

(b) Use part (i) to find two linearly independent solutions of this system. 


(c) Write down the general solution of the system. (You may assume that it is a linear 
combination of the two solutions you have found.) 


(ili) Can you find a solution of the form specified in (i) for the system 
X, =X 
X= —2x, + 3x,? 
(iv) Solve the system 
X, = Xz 
X2 = —4x, + 4x, 4+ @, 
using the steps outlined below. 
(a)In part (11) you calculated an eigenvector b of the matrix 


2 a 


and a vector a such that (B — Al)a = b. Form a matrix P whose columns are b and 
a (in that order). Show that 


PP 8F = ; 
0 2 


(b) Make the substitution x = Py in the given system, to obtain 
yi =2y, + y2 
yo = dy. + e. 
Solve the second of these equations for y2, and hence solve the first equation. 
(c) Hence solve the given system. 
[Solution on p. 40] 
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Appendix 1: Solutions to the exercises 


Solutions to the exercises in Section 1 


1. (i) Differentiating the components of 
g(t)= [2t7 4 ~~ e*]? we obtain 
a =[4t 0  26%]7. 
(ii) Let a= [a, a> a, |’. Then 
See =fae 2c >. age Y. 
Thus 
. t= [ae ee oe eT. 
= Ae“ [a, ene re 


= Age”, 
(This is just the result we would have hoped for.) 


2. Systems (1), (2) and (3) are linear, (4) is not. None of (1), 
(2) and (3) are homogeneous. (1) and (2) are constant- 
coefficient, (3) is not. (1) and (3) are first-order; (2) is second- 
order. 


3. (i) It is helpful to rewrite the system as 
X, — 3x, = —2x, —x) 
x3 = X1 +X2+%X; 
—X,+ X,-X3;= 0. 


In matrix form this is 


Ax = Bx (or Ax — Bx = 0) 
1 -—3 0 —2 -1l 0 
withA=]} 0 0 1 and B= i 1 1 
—1 1 -1 0 0 0 


(ii) It is constant-coefficient and homogeneous. It is first- 
order. 


4. (i) This system is A,x + Ax = h(t) where 


— ee ' 
a=|) a a=) 2 m= 


We need A; *. Since A, is 2 x 2, we can just write this down 
(see Unit 20, Subsection 4.1): 3 


oe ag Oe ee 
A;'=-43 = 
: {a F -4 
In normal form the system is 


x = —Ajy!A>x + Ay h(t), 


where 
her 27115 0; tts 
Se ee 10) tee 


1} 4 
A;'h(t) = : e | 
eee 


and 


So the system, in normal form is 


c —1 me st + % 
ct gage te 4 —4/ 


(ii) This system is A;x + Ax = h(t) where 


| 1 5 0 t 
m=ly ap Malt if wo=[o} 


The rows of A, are linearly dependent, so A, is singular and 
A; ' does not exist. Hence this system cannot be written in 
normal form. 


(iii) This system is not linear, and hence cannot even be 
written in matrix terms, let alone normal form. 


(iv) This is A,X = Aox where 


fo. y ee | 
ao - 
! I a ap oer F 1 


Now 
1 -1 4B + 
Avi = 4 = 
i a F ~4 
SO 
1} 1 1 
Aral at 
Sere ee ara 


In normal form the system is therefore 
x = A; 'Aox 
where A; ‘Ag is as above. 
5. The system is already in the normal form x = Bx where 
2 3 
= ; | 
We need to find the eigenvalues and eigenvectors of B. 


The eigenvalues are found by solving the characteristic 
equation det(B — AI) = 0. Now 


2-A 3 
det(B — 21) = 3 eat 
=(2—A)(1— 4) -6 
= 1? -3,-4 
=(A—4)(A + 1). 


So the eigenvalues are 4 and —1. 


To find the eigenvectors we solve 


Pa -allel-[0 
2 1—A]jv 0 
putting / equal to each eigenvalue in turn. 
Case 2 = 4 gives 
—2u + 3v=0 
2u — 3v=0. 


So an eigenvector corresponding to the eigenvalue 4 is 
[3 a, 


Case 4 = —1 gives 
3u + 3v=0 
2u + 2v = 0. 


So an eigenvector corresponding to the eigenvalue —1 is 
Geers oS 


We have two linearly independent eigenvectors, as required. 
Hence the general solution is 


3 1 
x= c,) > + c} ie 


xj = 3C,e" +Cre! 

x, = 2C\e" —Cse™ 
where C, and C, are arbitrary constants. 
6. (i) We have x = Bx with 


ate 
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Now 
S-—-A 4 
det(B — AI) = 
et ( Al) # 4 
=J?—5,+4 
= (A — 1) (A — 4). 


So the eigenvalues of B are 1 and 4. 
To find the corresponding eigenvectors we solve 
fe allel 
~i ={fiy 0 
putting 4 equal to each eigenvalue in turn. 


Case 4 = 1 gives u + v = 0 so an eigenvector corresponding 
to the eigenvalue 1 is [1 —1]?. 


Case A = 4 gives u + 4v = 0 so an eigenvector corresponding 
to the eigenvalue 4 is [—4 i 


There are sufficient eigenvectors, so the general solution is 
1 —4 
x = c,) ie + c} te 


x; = Cie — 4C,e* 
X= =—Cie + Coe™ 
(ii) We have x = Bx with 


Or 


5 -6 -6 
B= |-1 4 et 
3 -6 —4 
Now 
S5-A -6 —6 
det(B-—AI)=|-1 4-) z 
3 -6 -4-} 


= (5 — A) [(4 — A)(—4 — 4) + 12] 
+6[44+ 4-6] —6[6 — 12 + 32] 

—J3 +57 —- 81448 

“(*+ HOSP 


So the eigenvalues are 1 and 2 (repeated). 


To find the corresponding eigenvectors we solve 


5—A -6 —6 u 0 
—-1 4-) 2 vi={|0 
3 —-6 —4-Aliw 0 


putting A equal to each eigenvalue in turn. 


Case 4 = 1 gives 


4 -6 —-6llu 0 
—1 3 2{| v} =] O} 
3 -6 —S|iw 0 


Performing row operations we obtain 


0 6 2][ u 0 
-1 3 £2}ivi=lo 
a 5 0 
0 0 Ollu 0 
—1 0 Lijo}=101 
a oe 0 
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Therefore 
—u+w=0 
3v+w=0. 

So an eigenvector corresponding to the eigenvalue 1 is 
[-—3 1 —3)]!. 
Case A = 2 gives 

35 O46] fu] > Fe 

—1 P rate Bes gis! 
3 -6 —-6]|w 0 


This reduces to the single equation 
=—u + 2p + 2w = 0. 


Hence there are two linearly independent eigenvectors, for 
example [2 0 1}* and (2 1 ue 


We have sufficient linearly independent eigenvectors, so the 
general solution is 


=3 2 2 
x= ar Ljeo+ C, 0 je” = C3 ] ae 
=a 1 0 


or 
xX; = —3C,ée + 206" i 2C3e" 
X2 = Cie + Ce 
Xx3= —3C,e = Cre. 


(iii) We have x = Bx where 


oe eee 
det(B — il) =| | 


—1 —A 
=174+2442. 
So there is a conjugate pair of eigenvalues —1 + i. 


The equation for the eigenvector corresponding to the 
eigenvalue —1 + iis 


Po 4 


(—1—i)u+2v=0 
—u+ (1 —i)jv=0. 


The first equation is (1 + i) times the second and so both 
give 2v = (1 + i)u. Thus an eigenvector corresponding to the 


or 


eigenvalue A= —1+iisa=([2 1+ £i]". (The eigenvector 
corresponding to the eigenvalue 1 = —1 — i is therefore 
a= (2 1 — i]’, but this need not concern us here.) 


Now 


ac =| 2 Jes 
1+i 


= 2(cos t + isin t) 
eB (1 + i) (cost + isin t) 


a 2cos t 2sint 
=e ; = : ‘ 
cost — sint cost + sint 


Re(ae“’) = e™! or ; 
cost — sint 


Im(ae“) = e™' gases 
cost + sint 


Thus, by Procedure 1.2(a), we can write down the general 
solution in the form 


| 2 cost 2sint 
x=e'' C, : + C, : 
cost — sint cost + sint 


that is ~ 
x, =e '(2C, cost + 2C, sin t) 
Xz =e '((C, + Cz)cost + (C2 — C,)sint). 


Solutions to the exercises in Section 2 


1. Note that a good deal of the calculation has already 
been done in Example 1. We have x = Bx + h(t) where 


2 3 
B = 5 1 as in the Example, but now 
Nit) = fe ee. 


Hence Step 1 of Procedure 2.1 is already done. We have 


os 2 4-0) 
P=|) a . BP =| Py 
ae 
pa(t #] 
. Peers. | 


We next put x = Py, to obtain 
y =P 'BPy + P“'h(t), 


that is 
Pcl lb Galle] 
y2 0 —1{))2 $ —#|| 4e” 
or 
y, = 4y, + e* 


hf Sete? or air 2e”". 


Each of these equations can be solved using the integrating 
factor method. Writing the first equation in the form 


y—-4n= a? 


we see that the integrating factor is e “*. Multiplying both 
sides by this gives 


—4t a we Ga —2t 
—(e i=¢ “ee 
dt e 
Xe) ey = forma = —4e°%4C,, 
1.€. y, = ee nae te! 


Writing the second equation in the form 
y2 = nae Sea —2e, 


we see that the integrating factor is e’. Multiplying both sides 
by this gives 

d 

—(e'y,) = —2e*, 

7h! y2) é 


SO 
ey. = —$e* + C2, 

i.e. 

t 


y2 = C,e~' — $e". 


The solution of the original system is x = Py; that is, 
Se ae ie 1 || C,e** — 4e?! 
Sof te whee ae 


x; = 3C,e* + Ce * = 13! 


x2>= EL BOP ag = Ce" =— 47", 


or 


MST204 22 Solutions 


2. We have 
x = Bx + h(t) 


where B = = . and h(t)= [48 91e']7. 


Solving det(B — AI) = 0 gives the eigenvalues of B. Now 


725° = {69 
det(B — AI) = 
ee ee er 
= (2 + A)? — 10? 


= (A+ 12)(A — 8). 
So the eigenvalues are — 12 and 8. 


Eigenvectors corresponding to the eigenvalue 8 satisfy. 


—10 10 }juj} |0 
10 -10}{v} [of 
that is u = v. One such eigenvector is [1 if: 


Eigenvectors corresponding to the eigenvalue — 12 satisfy 


fo wollel-[o} 


that is u + v = 0. One such eigenvector is [1 —1]?. 
1 1 
Let? = | Then 
1 -1 
—-1 -1l 5 5 
Pola aj = 
= i k 4 


and 


et 
BP = 
. F of 


| Putting x = Py in the original system gives 


y =P 'BPy + P“'h(t). That is 
yi| | 8 0 |) 
Pl-Lo -alb}+L 


Yr. = 8y, +244 Ye 
y2 = 12y, + 24 — Bet. 


We can solve these equations by the integrating factor 
method. We obtain (omitting the details) 


y, = Cye™ — (3 + He’) 
y2 = C,e-'* + (2 — Fe’). 
Then 


sel! 1|[ C,e* — (34+ He’) 
ee ee 2 os eS eee 


| 


I[sie} 


Ni bo 
Sn) Ley 


or 


or 
x, =Cje + 0c 1 = ee 
w= Ce — Ce 4 5 = 3: 

3. (i) As in Example 2, we can write down recurrence 


relations corresponding to Euler’s method directly. We 
obtain 


Xirti1 = X1, + h(3rhX2, + 4) 
Xo r+ aes X2, + h(rhX,, Bre X2, a e”"). 


These are already in a form suitable for generating solutions, 
although it is a little neater to rewrite the second equation: 


Xo r+ = ae. Oe 5 (1 — h)X>, baa he”™. 
We have also that X; 9 = 5, X29 =2. 


(ii) With h = 0.1 we need to calculate X,, and X, > to 
approximate x, (0.2) and x,(0.2). Putting r = 0 in the 
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recurrence, we obtain 
X1,=5+01 x4=5.4 
X2,=0+09 x 2-—0.1 = 1.7. 
Putting r = 1, we obtain 
X12 = X,,4+01(00.3X,, 4+ 4) 
=54+0.1 x 4.51 
ae e 


and 
X22 = 0.01 x 5.44+ 0.9 x 1.7 —0.1e%! 
== 1.473. 


Thus we obtain the approximations x, (0.2) ~ 5.851 and 
x2(0.2) ~ 1.473. 


Solutions to the exercises in Section 3 


1. To calculate the eigenvalues of A, we first form its 
characteristic equation: 


»_\s-4 8 
$  -§-3 


= (A +3) — 
= (A+ 1)(A + 3), 
giving the eigenvalues —1, —4. 


To find the corresponding eigenvectors we solve 


—3-A § u| [0 
3 —j-—Allv Oo; 
putting J equal to each eigenvalue in turn. 
Case 4 = —1 gives the equations 


gu + $v = 0 (twice). 
Thus u= —vandso [1 —1]? is an eigenvector 
corresponding to the eigenvalue — 1. 
Case 1 = —% gives the equations 

—gu + gv =0 

gu — $v = 0. 

Thus u = v and so [1 1]’ is an eigenvector 
corresponding to the eigenvalue —4. 


2. The pair of differential equations can be written Ff = Ar 
where 


4 
ant i} 


Let us first find the eigenvalues and eigenvectors of A. To 
find the eigenvalues, form the characteristic equation. 
6 


eg 
0= 1A — at|=| : 
3 


§ 

oe og 
=(A+P)A+ 4) -% 
= A? + PA + ep 
=17+514+4 
=(A4+4)(A + 1). 

Thus the eigenvalues of A are —1 and —4. 


The eigenvectors corresponding to —1 are given by 


—42 4 4 8 Hy 40 
2 —,iiio| for 
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Thus v = 3u and so an eigenvector corresponding to the 
eigenvalue —1 is [1 31", 


The eigenvectors for —4 are given by 


3" sale) -Lel 
7 —P+4][v] Loy 
That is 
gu+$v=0 
Su+48v=0. 
Thus u = —2v.and so an eigenvector corresponding to the 
eigenvalue —4is[—2 _1]?. 


Since the eigenvalues are negative, the object will describe 
simple harmonic motion along the corresonding eigenvectors. 
The frequencies will be the square roots of minus the 
eigenvalues. So there are two possibilities: 


(i) motion in the direction of the vector [1 3]7, with 
angular frequency 1; 
(ii) motion in the direction of the vector [—2 1]7, with 


angular frequency 2. 
3. (i) We have X = Bx where 
—— k it 
1 0 
From the characteristic equation 
—A 1 
1 -A 


The eigenvalues of B are 1 and —1. 


0= =j? -1, 


The eigenvectors corresponding to the eigenvalue 1 satisfy 
—1 L}ju} |0 
1 —-1]}jv] [oy 


So one such eigenvector is [1 Lr. 
The eigenvectors corresponding to the eigenvalue —1 satisfy 
alle] [o] 
1 1 |i ov 0; 
So one such eigenvector is [1 —1]’. 


By Procedure 3.2 the general solution is 


1 
i= [|e + D,e"') + |: cost + D, sin t) 


or 
x; = Cie + Dye“' + C, cost + D,sint 
X2 = Cie + Die‘ — C, cost — D, sint. 


(ii) We must first divide through the second equation by 2, 
to obtain X = Bx with 


—-10 -6 
B-| —6 I 


From the characteristic equation 


a (1 eke OR 
es —6 1) = 
= (A + 10)? -- 6? 


= (A+ 16)(A+ 4) 
the eigenvalues of B are —4 and — 16. 
The eigenvectors corresponding to the eigenvalue —4 satisfy 
—6 -—6|]|] u 0 
= a : Bi 


So one such eigenvector is [1 —1]?. 


The eigenvectors corresponding to the eigenvalue — 16 satisfy 


& sifu) fo 
—6 6j}iv}| [Of 
So one such eigenvector is [1 i)’. 


By Procedure 3.2 the general solution is 


1 
x= ‘lc: cos 2t + D, sin 2t) 


1 
5 He cos 4t + D, sin 4t), 


or 


x, = C, cos 2t + D, sin2t + C,cos 4t + D, sin 4t 
Xx, = —C, cos 2t — D, sin2t + C,cos 4t + D, sin 4t. 


Solutions to the exercises in Section 4 


1. In each case we must write the equation in the form 
A,X + Ax + A3x = 0, and see whether the matrices A,, A, 
and A; Satisfy the specified criteria. 


(i) Here 


ree 2 -1 
es es re 3 
Each of these matrices is symmetric. A; has eigenvalues 1 


and 3. We must find the eigenvalues of A, and A3. 


For A, we have 


det (A, -ay=P >? ae 
= (3 — A)? — 2? 
= (A —1){(A — 5). 


So A, has eigenvalues 1 and 5. 


For A; we have 


det(A, — Al) = es an ; 
= (2-4)? —(-1)/ 
= (A — 1)(A — 3). 


So A; has eigenvalues 1 and 3. 

Hence all the eigenvalues of each matrix are positive. 

So the conditions of Theorem 2 are satisfied in this case. 
(ii) Here 


oe og 
A, = rn 
eG} St 
22H) 
A; = 
lol 


Each of these matrices is symmetric, and A, and A, have 
positive eigenvalues. However A; has 0 as an eigenvalue, so 
the conditions of Theorem 2 are not satisfied in this case. 


[In fact, one can see that x, = 0, x. = a (where a is a non- 
zero constant) is a solution of this system—and is not a 
transient. | 


(iii) To compare with the form in Theorem 2 we need to 
bring all the terms to the left-hand side of the equation . So 
here 


—2 int 
A, — IL, A> = I, A; -| 1 | 
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Each of these matrices is symmetric, and A, and A, have 
positive eigenvalues. However, _ 


det(A, — al) = (—2.—4y,- C1 
=(-(- 3-2 


Hence the eigenvalues of A; are negative (—1 and —3), and 
the conditions of Theorem 2 are not satisfied in this case. 


2. Ifx, = —1, x, = —t, then Xx, = X, = 0, and the equations 
are satisfied. 


From Solution 3(i) in Section 3, the general solution of 
X, =X, 
X2 =X 
—the associated homogeneous system—is - 
x, = Cie + Die '+ Ci, cost + Dz sint 
Xx, = C,e + D,e ‘ — C, cost — D,sint. 
So the general solution of the given system is 
x, = Cie + Die '+ C, cost + D,sint,— 1 
x, = C,e' + D,e'—C,cost — D2, sint — t. 
3. (i) We have a solution of the form x = Re(ze?") if 
—4z, + 4z, + 2z, =6 
—47, +2, + 922 = —2i. 
That is, if 
22526 
Z, + 5z2 = —2i. 


Thus z, = 3, z; = —2i — 15 and so the required particular 
solution is 


x, = 2sin2t — 15cos2t, 
X2 = 3cos2t. 


(ii) We have a solution x = Re(ze*) if 


—2Z, — 22, + iz; + 2; — 322 = —i 

—3z, — 22 + 2iz, + 22; + 22 = 14 21. 
That is, if 

iz; — 522 of 


—z, + 2iz, =1+2i. 


Adding i times the second equation to the first gives 
—7z, = —2,s0 


ae 
Z2 = 7. 


From the second equation 


4i 
z= — (1 + 2i) 


The required particular solution is therefore 
x; = —cost + 4sint, 
Xx. = ¥cost. 


Solutions to the exercises in Section 5 


1.4) The system can be written x = Bx with 


ee 
B = 
& es 


To find the eigenvalues of B we solve the characteristic 
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equation det(B — AI) = 0. We have 


Pi 5 
det(B — 21) = = 
=j2_1-6 
= (A — 3)(A + 2), 


so the eigenvalues of B are 3 and —2. 


The eigenvectors corresponding to the eigenvalue 3 satisfy 
5 —Silu 0 

Hs - i 4 < B 
Thus u = v and so an eigenvector is [1 ET. 
The eigenvectors corresponding to the eigenvalue —2 satisfy 

10 —S]iu 0 

i i 2 Bi 

Thus v = 2u and so an eigenvector is [1  2]*. 


By Procedure 1.2 the general solution of x = Bx is 


1 1 
x= cal te + cal fe 


That is 
xX, = Spe + Ce 
x2 = Ce + fee; 


(ii) This can be written x = Bx + h(t) where B is the same 
as in (i) and h(t) = [2e' —e*']’. The general solution can 
be found using Procedure 2.1. 


The eigenvalues and eigenvectors of B are given in (i). 


We form the matrix 


and put x = Py. This gives 


2e 
y =P 'BPy + | 7 


3 0 2 -1 2e' 
: i a, u * Leh 
That is 
¥, = 3y, + 4e' +e” 
V2 = —2y2 — 2e — e*. 


Using the integrating factor method we find that the- general 
solutions of these equations are 


y= Ca ae e2! 


y2= rs = ge! oom te7" 


So the general solution of the original system is x = Py, that 
is 


x= Cie" + Coe" = $6 — 22" 
xo Cie" + 2C.e" * — te — de" 


(iii) This can be written x = Bx where 


oS cae | 
10 —6 


The eigenvalues and eigenvectors of this matrix are complex. 
One eigenvalue is A = 2i, with corresponding eigenvector 
a= [4 6 — 2i]'. 


Now 


Pp 4 - 4(cos 2t + isin 2t) 
a = e = F 5) 
goer oe” (6 — 2i) (cos 2t + isin 2t) 
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SO 
4 cos 2t 
R At) _ 
= hi cos 2t + 2sin m 
4sin 2t 
At 
oe F sin 2t — 2cos sf 


By Procedure 1.2(a) the general solution is 

x = C,Re(ae”) + C,Im(ae“’), 
or 

X,; = 4C, cos 2t + 4C, sin 2t 

X2 = C,(6cos 2t + 2sin 2t) + C,(6sin 2t — 2cos 21). 
(iv) In matrix form this can be written ¥ = Bx where B is 


the same matrix as in (i). The eigenvalues and eigenvectors of 
B are given in (i). By Procedure 3.2 the general solution is 


: 1 
x = [cre + D,e~¥3") 


+ | |e: sin ,/2t + D2 cos ,/2t) 
or 
xX,= C,e%3! i Die v3" + C, sin pt + D,cos we s2 
X= C,ev3" + Dye V¥3t + 2C, sin 2 + 2D, cos «ft: 


(v) The general solution is the sum of the complementary 
function, which we found in solution (iv), and any one 
particular solution. We can find a sinusoidal particular 
solution by applying Procedure 4.2. 


Putting x = Re(ze?") leads to the equations 
—4z, = 8z, —5z, -i 
—4z, = 10z, — 7z, + 2, 

which have the solution 
= -F—- Yai 22 = — - Fi. 


Hence a particular solution of the system is 


xX, = —Fcos2t + sin 2t 
x2 = —+7cos 2t + ssin 2t. 


The general solution of the system is therefore 
xX, = C,ev3! = D,e7 v3" a 5 C, sin /2t + D, COS wee 
— Fcos 2t + sin 2t 
%5 == C,ev3! + Die v3! + 2C, sin sat + 2D, cos at 
— cos 2t + 3 sin 2t. 


(vi) This equation can be written x = Bx where 


=! 2 0 
B=j|-1 -3 0}. 
0 1 —4 


To find the eigenvalues of B we solve the characteristic 
equation det(B — AI) = 0. Now 


cae 2 0 
gett he i 0 
0 1 se er 
mee area 2 
es 
(an ap 


= (—4 —j)(A* +. 444 5). 
So the eigenvalues of B are —4 and —2 +i. 
The eigenvectors corresponding to the eigenvalue —4 satisfy 
3 2 O}} u 0 
—1 1 0 
0 1 0 || w 0 


So u=v=0,and aneigenvector is [0 0O 1)’. 


The eigenvectors corresponding to the eigenvalue —2 + i 
satisfy 


1-i 2 0 u 0 
—-1 -l-i 0 v| =| OL 
0 1 —2—-—i]|w 0 
That is 
(1 —i)u+2v=0 
—u—(1+iv=0 
v—(2+i)w=0.- 
The first equation is (i — 1) times the second equation. The 
second equation gives u = —(1 + i)v and the third equation 
gives w = v/(2 + i) = 4(2 — i)v. Thus an eigenvector 
corresponding to the eigenvalue 4 = —2 + i is 
a= [-(1 +i) 1 4(2 — i)]’. 


The general solution of the given system can be found using 
Procedure 1.2(a). We have 


—il-i 
ae*=e * 1 |(cost + isint) 
3(2 — i) 
—cost + sint —sint — cost 
=e 7! cos t +i sin t é 


4(2cost + sin t) 4(2 sin t — cos t) 
Se) 
—cost + sint 
Re(ae”’) = e 7! cost 


3(2cost + sint 


—sint — cost 


Im(ae“) = e 7! sin t 


+(2sin t — cos ft) 


The general solution is therefore 


0 —cost + sint 
x=C,|Ole*+4+C, cost em 
1 4(2cost + sin t) 
—sint — cost 
+C; sin t ” Sad 


+(2sint — cost) 


or 
x, = C,e° **(—cost + sint) + Ce *(—sint — cost) 
x, = C,e “cost + C,e “sint 
XxX, =C,e “ + 4C,e— “(2 cost + sint) 
+ 4C3e 7'(2sint — cost). 


(vii) We first write the equations in standard form. To do 
this, we separate the derivative and non-derivative terms: 


Xi +X» = 4x, + 4x, 
X = 3x oe Mg = ey 
X3 = 4x3. 
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and then write this in the matrix form Ax = Cx where 


1 1 0 4 4 0 
A=|0 1 0 and C43 1 1 
0 0 1 0 0 4 


The normal form of the equation is 
x=A “Cx 


The inverse A~’ can be found using row operations (see Unit 
20, Subsection 4.2). We have 


A I 
1 1 0 1 0 O|R, 
0 1 0 0 1 0|R, 
0) Co 0 0 1jR; 
R, —R,/ 1 0 0 1 -1 0 
0 1 0 0 1 0 
W) 0 1 0 0 1 
I A= 
Thus 
1 -1 
A-'=/0 1 0 
0 0 1 
So the normal form is 
1 -1 0||4 4 0 
x= |0 1 0||3 1 1 |x 
0 0 1||0 0 4 
1 3 -1 
==} 3 1 13% 
0 0 4 


We next need the eigenvalues and eigenvectors of 


1 ete 
B= |3 1 ae 
0 0 4 


These were calculated in Example 1 of Section 3 where we 
found 


eigenvalues corresponding eigenvectors 


By Procedure 1.2 the general solution is 


1 1 0 
x =C,e ~* |—-1 |+ Cre“|1| + C3e*'| 1 |, 
0 0 3 


Or 
% 206°" 26,2" 
— C +24 eal 41 
X27 = —C€3\e a Cre + C3e 
X3= fae a 
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2. The system can be written in the form 
A,X + Ax + A3x = h(t) 
with 


Each of the matrices A,, Ay and A; is symmetric. The 
eigenvalues of A, are 3 and 5 for 

(det(A, — AI) = (4— A)? — 17 = (A -—3)(A— 5), 
and those of A, are 3 and 3— all positive. Those of A; are a 
little harder to find. We have 


det(A; — Al) = (3 — A) (6 — A) — 16 
= 47 —9) 4+ 2. 


So the eigenvalues of A; are 4(9 + ./97 — 8). Written in this 
form we can see that these eigenvalues are positive, since 


oo FD: 


Hence we may deduce from Theorem 2 of Section 4 that the 
complementary function of the given system is a transient. 


There is a unique sinusoidal particular solution of the given 
system. Every solution of the given system converges to this 
same solution as t becomes large. Let us now calculate this 
solution. 


39 


To do this we put x = Re(ze*"). This yields the equations 
—16z, — 4z, + 6iz, + 3z, + 4z, =1 
—4z, — 16z, + 61Z2 + 47, + 6z, = —1. 


That is 
(6i — 13)z, = 1 
So 
1 6i + 13 
41; >= ee 
6i — 13 205 
soe i(6i+ 10) —6+ 10i 
ee 
The steady-state solution is therefore 
X, = —395(13 cos 2t — 6sin 21) 
X2 = —746(6cos 2t + 10sin 2r). 


3. The relevant recurrence relations, for step-length h, may 
be written down directly: 


Xirtt =X, +hX2, 

Xort1 = X22, +hX3, 

X3r41 = X3,+h(X,, + 3X2, — rhX3, —cosrh) 
Where A, 4 = 2, -X39 = —1; X35 =4 | 


(Here X, = [X,, Xo, 
value of x(rh).) 


X3,]' approximates the true 


40 


MST204 22 Solutions 


Appendix 2: Solutions to the problems 


Solutions to the problems in Section 5 


1. (i) We have y=x,,X, =X2, X= x3, and 


X3 =X, + 3x2 — tx3 — cost. (1) 


We can express the X3, X;, X2 and x3 in (1) in terms of y as 
follows 


Xi =y), 
X2=X,=y, 
X3 =X2=), 
X3=)y. 
Substituting these into Equation (1) gives 
y=y+t 3y — ty — cost, 
or 
y+.ty — 3y — y = —cost, 
as required. 


(ii) We are given that y satisfies ‘jy’ + 4ty = sin 2t. To reverse 
the process in (i), we introduce some new variables. We 
define x,, x, and x3 by the equations 


Xi =y 
X2=y 
X3;= j. 


Then we have © 
X; =X 
and 
- By =X. 
Also x3 = y; so from the given differential equation, 
X3 = sin 2t — 4tx,. 


Thus the given equation is equivalent to the first-order 
system: 


X, = X2 
X2 = X3 
X3 = sin 2t — 4tx,. 
The initial conditions give 
X:(0)=1, x2,(0)=2, x3(0)=0. 


We can use the method described in Subsection 2.2 to write 
down recurrence relations to generate approximate solutions 
to this system: 


Xiyti = X1,+hX2, 
X21 =X,,+hX;3, 
X3r41 = X3, + A(sin 2rh — 4rhX, ,). 
with X10 =1, X29 =2, X39 = 0. 
If we use these recurrence relations, then X,, will provide an 


approximation to the solution y(rh) of the original 
differential equation. 


(The method described in this problem generalises the 
method of numerical approximation for a single second-order 
differential equation described in Unit 6.) 


(iii) If we putx,;=x, x,=y, x3;=X, X4=Yy, we have 


X1 = X3 
X2=Xq 
—kx, 
X3(=X)= — A(x3 + x4)'/?x3 
xt + x3 
—kx, 
. <i ee a 2 2 1/2 
X4(=y) = = ax +X) Xe. 
Xi + xs 


The following recurrence relations are obtained if we apply 
Euler’s method to each equation in turn (with step length h). 
Xir+t = AX» + hX3, 
Xo r+ = Xo, i hX 4, 
kX,, 
Kaper = Xa, — A = 


ats ee M(X3,, + X4,)!?Xs,) 
3 BF Ze 


kX>, 
Xarti =Xay — | ete + A(X3, + X2,)!?Xe, | 
2. 8) 
x = ae” + bte*! (1) 
then 


x = jae” + b(e* + Ate). 
So (1) is a solution of x = Bx so long as 
(Aa + b)e* + bAte* = B(ae* + bte“). 
This is the case so long as 
da + b = Ba (equating coefficients of e*“), (2) 
and Ab = Bb 


That is, (a) b is an eigenvector of B, with eigenvalue / (from 
(3)), and (b) a satisfies (B — AI)a = b (from (2)), as required. 


(equating coefficients of te“’). (3) 


(ii) (a) We have x = Bx with 


[2] 


The eigenvalues satisfy 
—A 1 
—4 4-—; 


Thus 2 is the only eigenvalue. The corresponding 
eigenvectors satisfy 


—2 L}ju} {0 
—-4 2j}iv} [oO] 
Thus v = 2u, and so every eigenvector is of the form 


kK 2kjF=kf1 27°. 


To be able to apply Procedure 1.2 of Section 1 we require 
two linearly independent eigenvectors, but this is impossible 
since all eigenvectors are multiples of [1 fee 


c= 


=4?-4144+4= (4-2). 


(b) Since [1 2]’ is an eigenvector corresponding to the 
eigenvalue 2 we know that 


1 
ae Ba 


is one solution of x = Bx. We can use part (i) of the problem 
to find a second linearly independent solution. To do this we 
need to find a vector a satisfying (B — AI)a = b, where b is 


the eigenvector, [1 2]". That is, if a = [a, a] we 
need \ 
—2 a, ee l 
= Geen 5 Me eee 
That is 
—2a, + ay — 1 
—4a, + 2a, = 2, 


These two equations are equivalent, so there is a wide choice 
of solutions. One is a, = 0, a, =1 giving a = [0 Be 
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With these values for a and b, we know from part (i) that 


0 | 
X= fe + LS fe" 


is a solution of x = Bx. 


(c) The general solution is 
1 0 1 
x = c,) fe + c,(| Je" + fe 


=e Ce" Oe 
= 10 eS Ce" 20 te 


(To get this we assume the fact that the general solution is a 
linear combination of the two solutions we have found.) 


Or 


(ili) We now have x = Ax, with 


eos 0 1 
ape aa 
The eigenvalues of A satisfy 


1 


Thus the eigenvalues are 2 and 1. 


The eigenvectors corresponding to the eigenvalue | satisfy 
= aftel-Lo] 
—2 2]|lv 0} 
Thus u = v, and so an eigenvector is [1 es 


Let us try to find a solution of the form specified in part (1). 
We would need a vectora=f[a, a>]" satisfying 


2 allel-Lih 


That is 
a,—a,=1 
2a eS 2a, =a f. 


These equations are inconsistent and have no solutions. 


If we proceed similarly for the eigenvalue 2 we again arrive 
at a set of inconsistent equations for a. So in this case no 
solution of the form specified in part (i) can be found. This is 
just as well, as we have enough independent solutions 
already! 


(iv) (a) In part (ii) we had b= [1 2\',a= {0 1]’. 
[You may have different vectors here—there was more than 
one valid choice of these vectors. ] Thus 


41 
Then 7 
es 1 0 
2 iy 
and so 
P-'BP — 0 0 ri 0 
eae id ano 4 


(b) The given system can be written 


|e} 


Putting x = Py gives 


0 
Py = BPy +| | 
e 


x = Bx 


ob 


SO | 
. -1 _1| 9 
y a P BPy “+ P ot 
1.€. 
= 2 | ie 0 
— 0 2 y e 
Thus 
V1 =2y, + yr 
y2=2y.,+ e. 


We can solve the second equation, by the integrating factor 
method to obtain 


yz = Cye* — et 
Substituting this into the first equation gives 
V4 =—— 2y; + C2 aot e'. 


We can now solve this by the integrating factor method, to 
obtain 


yy, = Ce + Cite + e. 


(c) Now x = Py, so we have 


1 0 }{ Cre” + C, te + 
x= 
2 1 Cie" —¢ 


Thus 
xX, = OY a Cte" +e 
Xo = 2C,e" + Cye"(1 4+ 21) + e 


is the general solution of the given system. 
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